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Разработана математическая модель определения характеристик многомодового полупроводникового лазера с шероховатой 
поверхностью активного слоя. Показано, что указанные шероховатости вызывают радиационные потери волноводных мод, а 
также приводят к формированию стохастических межмодовых связей внутри модовою спектра. В частности, установлено, что 
данные эффекты оказывают на параметры лазера воздействие, аналогичное увеличению коэффициента Петермана. Определе- 
но количественное соотношение между уровнем спонтанной эмиссии и статистическими характеристиками шероховатостей оп- 
тического волновода. 


Постановка задачи 

Активные оптоэлектронные волноводные 
структуры: полупроводниковые лазеры (ППЛ) и 
усилители, а также допированные ионами эрбия 
волоконные световоды широко применяют в опти- 
ческих системах связи, системах хранения и записи 
информации, технике высокоточных измерений и 
других областях современной техники. Производи- 
тельность данных систем определяется свойствами 
используемых в них волоконных и полупроводни- 
ковых лазеров, и, прежде всего, их шумовыми и 
модуляционными характеристиками. На указан- 
ные параметры заметное влияние могут оказывать 
технологические погрешности производства ак- 
тивных структур. Так, в §аіп цшйеА РР-лазерах 
(Фабри-Перо лазеры с рефрактивным волноводом) 
резонатор формируется в виде отрезка одномодо- 
вого оптического волновода (О В), представляюще- 
го собой тонкий активный слой двойной полупро- 
водниковой гетероструктуры (рис. 1) [1-4]. 



Рис. 1 . Волноводная структура ППЛ: 7 - длина резонатора, 
ѵѵ - ширина полоски контакта 


Динамика спектрального состава излучения та- 
кого ППЛ, представленного набором продольных 
волноводных мод (ВМ), зависит от типа ОВ. В ли- 
тературе она описывается коэффициентом Петер- 
мана (Ре(етіапп-РакГог) К !р [1]. При этом считается, 
что как КШ ( геіаііѵе іпіетііу поізе), так и спектр из- 
лучения ППЛ полностью определяются пропор- 
циональным К !р уровнем спонтанной эмиссии, вы- 
зывающей стохастические флуктуации амплитуды 
и фазы вынужденного излучения. 

Однако, аналогичные флуктуации числа вы- 
нужденных фотонов 5 ЯІ в /-ой моде лазерного дио- 
да (ЛД) могут быть вызваны и другими механизма- 


ми. Подобное воздействие на оказывает, напри- 
мер, стохастическая связь между ВМ, обусловлен- 
ная их рассеянием на микронеоднородностях неи- 
деальной поверхности ОВ. В зависимости от техно- 
логии изготовления эта поверхность всегда харак- 
теризуется в той или иной степени выраженным 
случайным полем локальных микроскопических 
отклонений (шероховатостей) от среднего уровня. 
Указанные шероховатости могут оказывать влия- 
ние на параметры лазера. 

Прежде всего, шероховатости приводят к до- 
полнительным радиационным потерям энергии 
ВМ в резонаторе и соответствующему увеличению 
порогового тока лазера. Физический механизм 
указанных потерь связан с преобразованием ВМ в 
непрерывный континуум излучательных мод 
(ИМ). При этом межмодовая связь поддерживается 
набором элементарных гармонических решеток, 
образующих трехмерный энергетический спектр 
шероховатостей пленки ОВ (7(К). Соответствую- 
щий коэффициент радиационного затухания а г 
выражается через парциальные коэффициенты 
а(ѵ,р) как [5]: 

л-ѵі л 

а,. = | | а(ѵ,Р)с!ѵсір , 
о о 


ѵ, р - сагиттальный и азимутальный углы, опреде- 
ляющие направление ИМ; угол іЛ определяется 
числовой апертурой ( питегісаі арегіиге КА) ОВ: 
1ЧА=$іп(іЛ); 

а(о,р) = С(К(о,рШ(о,р)) 2 х 


х 
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Е л (х) - нормированный профиль ВМ; р(ѵ.(і) - ко- 
эффициент связи ИМ и ВМ. 

Другой механизм влияния шероховатостей ак- 
тивного слоя ЛД на его характеристики связан с 
формированием стохастических межмодовых свя- 
зей (8МС - хіосИазгіс тосіе соирІіп§) внутри медово- 
го спектра лазера. В силу стохастической природы 
этих связей, число фотонов рассеянных в /- ую 

из у- ой моды всегда 5-коррелированно с собствен- 
ным полем фотонов У, данной моды. Таким обра- 
зом, оказывается, что воздействие флуктуаций А8 ір 
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на 5; вполне аналогично воздействию на них фото- 
нов спонтанной эмиссии. 

Данный механизм межмодовой перекачки фо- 
тонов реализуется лишь при фазовом синхронизме 
взаимодействия, при котором волновые числа Д 
ВМ связаны с вектором решетки К из спектра ше- 
роховатостей 6'(К) соотношением: 

д-д+к=о. 

Учитывая различие частот ВМ, обусловленное 
особенностями самовозбуждения лазера, из по- 
следнего соотношения можно заключить, что меж- 
модовая перекачка возможна лишь в режиме гене- 
рации нескольких поперечных мод, когда частот- 
ный интервал А у, между ними кратен интервалу А ѵ, 
между аксиальными модами: 

Ду, = тАѵ, , 

где т - целое число. 

Последнее условие не может удовлетворяться в 
полупроводниковых ОРВ-лазерах. работающих в 
одномодовом режиме. В этих устройствах техноло- 
гические шероховатости волновода приводят лишь 
к дополнительным потерям и повышению порого- 
вого тока. Однако эти условия могут выполняться в 
полосковых РР-лазерах, в которых возможна гене- 
рация на нескольких поперечных модах. Резонан- 
сное усиление взаимного рассеяния энергии мод в 
данных условиях приводит к дополнительным шу- 
мам излучения как отдельных ВМ, так и Мойе Раг- 
Ііііоп ІЧоі$е, а также оказывает влияние на характер 
динамических искажений оптических сигналов. В 
литературе указанные эффекты остаются практи- 
чески не исследованными. 

Поэтому целью настоящего сообщения является 
представление соответствующей математической 
модели и результатов расчетного моделирования. 


Скоростные уравнения 


Для моделирования динамических характери- 
стик ППЛ будем использовать модель так назы- 
ваемых скоростных уравнений [1—4]. 

Рассмотрим активный /Ѵ-модовый волновод- 
ный слой полоскового ППЛ на основе двойной ге- 
тероструктуры с резонатором Фабри-Перо. Дина- 
мика модовою спектра в данной структуре описы- 
вается системой дифференциальных уравнений 
для модовой плотности фотонов А, /-ой моды и 
концентрации носителей заряда п: 


сіі 
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где К яі - коэффициент вынужденной рекомбина- 
ции, определяющий скорость создания фотонов 
/-ой моды за счет вынужденного излучения; К 5рі - 
коэффициент спонтанной рекомбинации, опреде- 
ляющий скорость создания спонтанных фотонов в 


/-ой моде; т рКІ - время жизни фотонов; К ш - пол- 
ный коэффициент увеличения спонтанного излу- 
чения; I- ток инжекции; е - элементарный заряд; 
V - активный объем ППЛ; г„ - время жизни носи- 
телей заряда. 

Прирост числа фотонов сі.% в /-ой моде за время 
ёі определяется усилением уже имеющихся, вкла- 
дом спонтанного излучения, а убыль - потерями в 
резонаторе за счет материального поглощения и 
вывода излучения за пределы ЛД. Последние про- 
цессы формально учитываются введением времени 
жизни фотона в резонаторе т рІІ /. 

77 " = Іг іа "‘ ~ 21 ® 

рк,і V 


где п ег - групповой показатель преломления; а т і - 
коэффициент затухания в резонаторной моде; К ь 
Я 2 - коэффициенты отражения торцов резонатора. 

В системе ( 1 ) параметр К яі определяется соотно- 
шением: 

,Г„ (3) 

где §і - функция усиления для /-ой моды; Ѵ еі - 
групповая скорость фотонов /-ой моды; Г, - коэф- 
фициент локализации /-ой моды. 

Скорость генерации спонтанного излучения 
выражается формулой: 

=и,Лм’ < 4) 


где п !р - коэффициент спонтанного излучения. 

Функция усиления § в общем случае включает в 
себя зависимость как от А, так и от и и задается сле- 
дующим образом: 
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где ЛИ - так называемое дифференциальное уси- 
ди 

ление; п А — пороговое значение концентрации ин- 
жектированных носителей; е - коэффициент сжа- 
тия усиления. 

Система уравнений (1) в рамках данной модели 
полностью описывает динамические характери- 
стики ЛД с идеальной поверхностью активного 
слоя. 


Вырождение лазерных мод и межмодовая связь 

Рассмотрим далее влияние микронеоднородно- 
стей (рис. 1) активного слоя ППЛ на режим генера- 
ции излучения. 

Активный слой ППЛ представляет собой уча- 
сток планарного волновода (рис. 1). Влияние нео- 
днородностей границы планарной волноводной 
структуры на её волноводные свойства достаточно 
подробно изучено в литературе [1-4, 6, 7]. Вслед- 
ствие этих шероховатостей лазерные моды обмени- 
ваются энергией в процессе распространения по 
ОВ. При этом фотоны, рассеянные в /- ую моду ока- 
зываются (5-коррелированными по фазе и поляри- 
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зации с собственными фотонами. Подобный меха- 
низм межмодовой связи аналогичен влиянию 
спонтанного излучения на параметры ППЛ и мо- 
жет являться причиной дополнительных шумов и 
уширения спектральной линии, ухудшения его мо- 
дуляционных характеристик. 

Для использования подхода [1, 2] для постро- 
ения модели ППЛ рассмотрим различия механиз- 
мов межмодовой связи в пассивном планарном 
волноводе и активном лазерном резонаторе. В пас- 
сивном волноводе все моды, как известно, имеют 
одну и ту же частоту. Это обстоятельство упрощает 
достижение межмодового синхронизма и является 
необходимым условием для формирования много- 
кратных межмодовых связей внутри модового 
спектра волновода. 

В лазерном же резонаторе условие самовозбуж- 
дения обеспечивается лишь для определенного на- 
бора частот: 



где п - показатель преломления заполняющего ре- 
зонатор вещества; а„ а у , а, - линейные размеры ре- 
зонатора; /, т, р — целые числа, характеризующие 
порядок резонаторной моды; с - скорость света. 

Отсюда видно, что все типы колебаний лазер- 
ного резонатора имеют различные частоты. Следо- 
вательно, в данных условиях не следует ожидать 
эффективного взаимодействия поля рассеянных 
фотонов некоторой /-ой моды с полями других мод. 

Однако из (6) следует, что при определенных 
значениях линейных размеров резонатора ППЛ и 
индексов /, т , р возможно так называемое вырож- 
дение мод, т. е. совпадение резонансных частот для 
поперечных мод разного порядка. Обычно, изме- 
нения линейных размеров порядка 3...5 % при из- 
готовлении ППЛ не контролируются [3, 4]. Следо- 
вательно, для некоторых экземпляров ППЛ можно 
ожидать резонансного проявления эффекта силь- 
ного взаимодействия лазерных мод и обусловлен- 
ного этим эффектом увеличения уровня шумов. 


Случай модового вырождения, при котором ча- 
стотный интервал между соседними поперечными 
модами 5ѵ Т кратен интервалу между продольными 
модами 5ѵ ь , показан на рис. 2. 

Рассмотрим условия возбуждения поперечных 
мод высокого порядка. Вблизи порога излучения 
лазеров поперечные моды высшего порядка обыч- 
но подавлены. Однако, при достаточно больших 
значениях мощности основной моды, в ряде кон- 
струкций ППЛ условие возбуждения поперечных 
мод низкого порядка выполняется вследствие про- 
странственного выжигания носителей заряда 
[1-4, 6, 7]. В соответствии с (6), эти моды генери- 
руются на частотах выше, чем частота основной 
моды. Рассмотрим, например, лазер с зарощенной 
гетероструктурой, работающий на длине волны 
А=1,3 мкм с шириной полоски \ѵ=3 мкм. Предпо- 
ложим, что основная мода и мода первого порядка 
ограничиваются шириной полоскового электрода. 
Приближенные выражения для пространственной 
структуры световых полей в резонаторе имеют вид: 

Е 0 = со5(К 0 2)$іп(Р2), Е 1 = 8 т(ЛТ 1 х) 8 Іп(/? 2 ), 

где К^л/\ѵ, К { &2 л/уѵ, ^=Мк/Ь. Пренебрегая далее 
дисперсионной зависимостью показателя прело- 
мления п, из волнового уравнения находим связь 
между модовыми индексами: 

(1 2 /'м 2 ) + {р 2 /ь 2 ) = (2я//с) 2 . (7) 

Учитывая малость члена (/ 2 /Ѵ) ур. (7) фактиче- 
ски определяет значение модового индекса р. Так 
для лазера с параметрами п= 3,5, А=1,3 мкм и 
1=250 мкм, находим: 

р и 2иТ/А = 1346. (8) 

Отсюда следует, что частотный интервал между 
продольными модами в рассматриваемом случае 
будет: 

А/;=с/2я1 = 0,171 ТГц. (9) 

Для указанного значения р частотный интервал 
А( т между поперечными модами определяется под- 
становкой в (7) значений модовых индексов основ- 
ной и поперечной моды первого порядка 1=1 и 1=2 
для тѵ=3 мкм: 



V 


ѵ 0 

Рис. 2. Вырождение лазерных мод 
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А/ т « Зс 2 /8иѴ/ = 1,327 ТГц. (10) 

В соответствии с (7), эта величина превышает 
межмодовый интервал продольных мод. Отсюда 
следует, что в рассматриваемых условиях основная 
мода и поперечная мода первого порядка с одним и 
тем же продольным индексом будут отделены друг 
от друга приблизительно восемью продольными 
модами. Для данной частоты разность в индексах р 
основной и поперечной моды первого порядка мо- 
жет быть определена из (7) как: 

Ар и ЪІР /2сѵ 2 р = 7,7. (11) 

Заметим, что при вариации Ар на 1/2, разность 
частот основной и поперечной моды изменяется 
ровно на половину продольного межмодового ин- 
тервала. Этому сдвигу соответствует изменение 
длины резонатора Ь, либо его ширины \ѵ всего на 
3 %. При изготовлении лазера столь малые измене- 
ния длины резонатора контролировать очень слож- 
но. В данных условиях при случайной вариации 
температуры, тока инжекции и т. д. возможно фор- 
мирование условий частотного вырождения и 
сильного модового взаимодействия лазерных мод. 

Полагая, что данные условия выполнены, оце- 
ним уровень межмодовой связи. При этом для опи- 
сания взаимодействия мод из области вырождения 
воспользуемся скоростными уравнениями ( 1 ) и ре- 
зультатами работ [5, 8]. 

Модифицированные скоростные уравнения 

В расчет будем принимать взаимодействие про- 
дольных и поперечных мод из рассмотренной вы- 
ше области вырождения. 

Связь лазерных мод в «холодном» резонаторе, 
обусловленная рассеянием на шероховатостях вол- 
новодного слоя, описывается уравнениями связан- 
ных мощностей [8]: 

ГІР ы 

= -2 а ,р, + х 4; с/( \а ; )[Р, -р,\ (12) 

]=иі*і 

где Р„ а і - мощность и коэффициент затухания і- ой 
лазерной моды; С(Акц) - энергетический спектр 
шероховатостей; <К - коэффициент связи между і- 
ой и у-ой модами, определяемый выражением: 

^ = | Аг (х) ■ Е і (х) ■ Е] ( х)(іх , 

где Е,(х) - амплитудный профиль /-ой моды, Ае(х) 
- флуктуация диэлектрической проницаемости. 

Учитывая далее, что групповая скорость Ѵ е ли- 
нейно связывает продольную координату і и время 
1 лазерной моды, и, представляя объемную плот- 
ность фотонов У через мощность, выразим динами- 
ку У через пространственную зависимость Р в (12): 

Р8_ = ]_РР_ 

Рі кѵ Рг 

где к - постоянная Планка. 

Преобразуем систему (1) с учетом вышеуказан- 
ных выражений к виду: 


мдг^-г^^). 


Рі еѴ т 




( 13 ) 


Полученная система нелинейных связанных 
дифференциальных уравнений в рамках созданной 
модели полностью описывает динамику модовой 
плотности и концентрации инжектированных но- 
сителей с учетом механизма межмодовой связи. 

Уровень шумов ППЛ количественно характери- 
зуется относительной мощностью шума (КЛѴ) [1,2]. 


кт = 2А /Ж (а,) 

(г ) 1 


(14) 


где А/ - полоса пропускания; \Ѵ Р (со) - спектраль- 
ная плотность мощности; со - круговая частота; 
(Р) — среднее значение мощности. 

Соотношение (14) выведено для некоторого 
обобщенного шумового процесса. При изучении 
конкретных механизмов возникновения шумов, 
например, спонтанного излучения, или шума мо- 
дового распределения ( Мойе Рагііііоп Шізё) необхо- 
димо учитывать специфику данных процессов. 

Следуя [1,2], оценку шумов ППЛ будем произ- 
водить на основе введения в систему скоростных 
уравнений (13) внешних шумовых источников, на- 
зываемых ланжевенскими источниками. Соответ- 
ствующие им флуктуации модовой плотности фо- 
тонов и концентрации инжектированных носите- 
лей обозначим, как //,(/) и Р„(() соответственно: 


РА, 

Рі 


+Р.Л*) 


Рп 

Рі 


/ гр 

еУ т „ 




(15) 


Свойства случайных процессов, описываемых 
ланжевенскими источниками, хорошо изучены 
[1,2]. Для проведения дальнейших расчетов будем 
считать, что силы Р{1) удовлетворяют условиям ста- 
ционарности и эргодичности, а их статистическое 
распределение подчиняется нормальному закону. 
Кроме того, их среднее значение будем считать ну- 
левым, а корреляционные свойства будем описы- 
вать 5-функцией. 

Флуктуации концентрации носителей, описы- 
ваемые членом Р„(1), обусловлены дробовым и ре- 
комбинационным шумами. Однако в [1] показано, 
что шумы ППЛ, в основном, определяются флу- 
ктуациями числа фотонов //,(?). На этом основании 
далее мы будем пренебрегать членом Р„(і). 

Кроме этого, так же как и в системе ( 1 3), ур. (15) 
дополним членом, отвечающим за межмодовую 
связь лазерных мод с вырожденными частотами. 

Для оценки КШ рассмотрим установившийся 
процесс генерации ЛД, при котором модовая плот- 
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гость фотонов и плотность инжектированных но- 
сителей п складываются из стационарного значения 
и шумовой составляющей с нулевым ср 

8,={8<)+А8„ 
п = (п) + Ап, 

(8,) и (п) в (16), (17) являются стационарными ре- 
шениями системы дифференциальных уравнений 
(15). 

Подставим в (15) выражения (16) и (17) и пере- 
ведем анализ в частотную область. Получим: 

О' со - 1 / т, )А8 І (со) = АР '' ,. (со), 


0-1/*. )Ли(ет) =-- у .Д 8 /со) 


і = і 


(18) 


со.. 


СО,, 


(16) 

ш $ р 

6> 

(17) 

со , = \\ 

Км 

Эл 


V 




Ѵт 

' 1 рЬ.і 

Далее следуя работе [1], получим окончатель- 
ные выражения для КШ основной моды и для шу- 
ма модового распределения (Мосіе РагШіоп /Ѵоме) 
соответственно: 

, ,2 т“І +со" 4К, ы п 

КІМ 0 =Н 2 (со) т^сі 2 , (19) 


с р „со,. 


(*5о 


где Л /,(<у) - Фурье-образы ланжевенских сил; т, - 
характеристическое время жизни /-ой моды, 
г/=//,+г/,.; т' - константа, определяемая соотно- 
шением: 


хм Р =- 


4 А/т 


Р* 


у (^/ 

'^І + ІШТ,) 2 ’ 


^ іщО.р \8 (] ! 

где у — эффективное время жизни фотона. 


( 20 ) 


7=1 


би 


где г, - время жизни носителей заряда. 

Среднее значение (А /,(&>)) определяется пре- 
образованием Фурье от автокорреляционной 
функции (18): 

(КіИ 2 )= 2 ^«,р(^)ѵ.г'- 


Значения А / ,(&>) и ДУ ( (й>) связаны между собой 
передаточной функцией Д(ю) [1]: 


А/ (ю) 

АС. ,(ю) 


ЯД®) = со,. 2 


І/т'+ую 

(./©) 2 +^ +1 

® 7 ., 


Выражения (19), (20) являются формальной ос- 
новой для моделирования шумовых характеристик 
лазерного диода и исследования влияния техноло- 
гических шероховатостей волноводного слоя на 
уровень КШ. 

На рис. 3 представлены результаты моделирова- 
ния динамики генерации излучения ІпОаА$ лазер- 
ного диода, полученные на основании системы ура- 
внений ( 13). На рисунке также показано установив- 
шееся значение 8,(1) фундаментальной и первой по- 
перечной моды. Пунктирные линии соответствуют 
динамической характеристике лазерного диода с 
шероховатой поверхностью активного слоя. Радиус 
корреляции и среднее значение величины шерохо- 
ватостей равно 1 нм, значение К !р равно 2. 



Рис. 3. Динамика излучения лазерного диода 


Сама передаточная функция для малосигналь- 
ного воздействия легко находится из (18) и опреде- 
ляется выражением [1]: 


Н,(со) = 


ДА. (со) 
АС,» 


1 


со. 


2 



+^-+1 

СОр, 


1/т '+ую’ 


где со ф со ЛІ — частоты релаксационных колебаний и 
демпфирования соответственно: 


Заключение 

Предложенная модель полупроводникового ла- 
зера с резонатором Фабри- Перо позволяет учиты- 
вать влияние технологических шероховатостей его 
активного слоя на характеристики излучения: уро- 
вень относительной интенсивности шума, порого- 
вого тока, спектральных характеристик и т. д. По- 
лученные данные показывают, что соответствую- 
щий механизм флуктуаций вынужденного излуче- 
ния лазера аналогичен механизму влияния поля 
спонтанных фотонов на шумовые и спектральные 
свойства лазерного излучения. Указанная связь мо- 
жет быть формализована в виде соотношений меж- 
ду приращениями коэффициента Петермана и ста- 
тистическими характеристиками шероховатостей 
поверхности. 

Установлено, что в лазерном диоде с микронео- 
днородностями активного слоя волноводные моды 
обмениваются энергией за счет их рассеяния на 
шероховатостях и стохастических межмодовых 
связей. Доля энергии, задействованной в процессе 
обмена, составляет около 30 %. 
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Показано, что в спектре стационарного излуче- 
ния лазера с идеальной поверхностью волноводно- 
го слоя полем паразитных мод можно пренебречь. 
Однако даже весьма малые шероховатости волно- 
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